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TR 166 - Vo! . IV 
I . INTRODUCTION 

The development of nozzles for hypersonic aircraft requires 
the ability to analyze the behavior of high temperature gases, 
often chemically reactive, flowing through ducts of compli- 
cated geometry. For this purpose, analyses such as those 
presented in the preceeding volumes of this report have been 
and will continue to be developed. However, the extremely 
complicated nature of these flow fields necessitates the a- 
doption of certain simplifying assumptions in the construction 
of the analytical models. It is therefore both necessary and 
desirable to carry out the parallel development of an experi- 
mental program aimed at assessing the accuracy and areas of 
applicability of the analytical results as well as at the 
acquisition of data for those physically interesting configura- 
tions for which suitable analytical tools may not yet be avail- 
able. 

Since the exact duplication of all gas conditions at the nozzle 
entrance station, especially the extremely high temperature 
levels, makes experiments both difficult and expensive, it is 
logical to investigate the possibility of simulating the actual 
nozzle flows with low temperature nonreactive gases. To pursue 
this goal, it is necessary to develop the similarity laws for 
nozzle flows and determine the parameters which must be dupli- 
cated. This effort has been carried out and the results of the 
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analysis are presented here. In addition to the development 
of these simi 1 ari ty . rul es , a number of cool gas flows have 
been considered to determine their suitability as substitutes 
for actual high temperature gases in an experimental program. 
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II . SIMILARITY REQUIREMENTS 

The proper simulation of a flow field requires that at any 
geometri cal ly similar point the Mach number and flow direction 
be the same as those which exist at the correspond!' ng point in 
the alternate flow and that the ratio of static pressure from 
one point to another in each of the flows be identical. The 
necessity to match the variation of p with 6 requires that the 
value of r used in the characteristics equations and defined 
for equilibrium chemistry as 

r = £( 9 P) 

P 3p S 

be duplicated at corresponding points. The requirement that 
the Mach number and flow deflection be identical at correspond- 
ing points requires that the initial conditions be similar in 
p and identical in M and 8. In addition, the variation of V/V^ 
ana a/ai with p must be the same, where and a-,- are the gas 
velocity ana speed of sound at the reference initial conditions. 

It should be noted that although the ratio V -j / a i must be the 
same for the nozzle and the model, the actual values V n - and a-j 
can be different since only the local Mach number must be the 


s ame . 
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The assumption is made that the chemistry is in equilibrium. 
Therefore, the conditions of similarity require that the varia- 
tion of pressure ratio (p/p-j), r and Cp with temperature ratio 
(T/T-j) be the same at geometrically similar points. This will 
assure that the variation of velocity ratio, Mach number and 
pressure ratio, as functions of temperature ratio are also the 
same. These conditions imply that the variation of three 
thermodynamic parameters with temperature must be closely simu- 
lated in the model, within the range of temperature ratio exist 
ing from the throat to the nozzle exit. 
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III . THERMODYNAMIC RELATIONS 

Consider a volume containing a unit weight of gas in thermo- 
dynamic equilibrium consisting of n moles and with a static 
enthalpy H per mole. In this case, the following thermodynamic 
relations are valid if the gas is a mixture of species i and n 1 - 
is the number of moles of each specie: 

In i = n (1) 

‘It’p ■ " c p ' T <tr> P < 2 > 

whereCp is the specific heat per mole at constant pressure. 

The flow is isentropic between discontinuities and is in chemi- 
cal equilibrium. Therefore, 

(|— ) p dT + (~)j dp = 0 (isentropic flow) (3) 

(|y)p dT + (ly)y dp = 0 (equilibrium chemistry) (4) 

The Gibbs free energy is given by 

G = H - TS ( 4 ) 

and 


dH 


T dS - Vdp 


(5) 
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because 


( *V) 


as: 

■ap- 


m)r 




but Equation (3) 


gives 


Therefore , 





The 


e Quation 


of state gi Ves 



Therefore 



( 6 ) 


O) 


( 8 ) 


( 9 ) 


( 10 ) 


( 11 ) 
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If w.e define 

wne rs .. 


_ /31n ni 

1 1 n VafFfV P r n jJ ' 


ana 


T fit? relatlda<‘.betw®eB' , ''p and .-p; i s -:-g ;» ?sa:'bp_.- 

3 In 


a T = (i-LD-E.) /^£.v£ 

s. " - * - • "7777° 

' v TT"? c ■ 


1 + Dj 


(13) 

( 1 3 .) ' 


ay must be the same at corresponding points or at the same 


9 C«T'Ue of temperature ratio (T ) . i n addition, since 

ow, 1 

R IT 

= £n ■{ H -j 


V 

p = 


nH 


w h >} t a .. W • i $ : t h e m o l e c « 3 a ? w e 4 § a t . a n d ■. n W • 

then 


(?0) 

(14) 


I Rr.addIU do-., , 


n c r 


{ l&i C 


if 

RoT 


, 3nH . 

( 3T P 

(IM,.. + 

a i 


^ ni G p-j + E 


n i H ,- 


/ d K •, D Vf 

v i. . * • n-r—r-;; 


°T. 

l 


(15) 

(f-i); 


where Dj^ is defined as 


a»< 


° T i = 

•V Lf 4 *-, fj 

/ O « -5 ... { & « \ 

* ...Afcrt* ft*.- 1 V---KS* » • 

4 -if- $■ ■ ’ 3 -I p 


31 n n is 

,ap, 

^ v .—/ T- P s. rr/ c 


o J 


\ i - ./ : 


Thus 

rVe® if wue- <lefn nt- 


°p: 


n i c Pi n -H ■ • 

£ — 1 + z n i n i 1 i 

n nT 

I c ; i n n\ 

\ « Tv.** ~ 7 >v • 


(17) 




c; v p- 


a:-’ 
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where 



• 


E n i 0 Ti = n D t 

• 

(18) 

The relation 

r V - 

between p and p is given by 




=- a 2 = aT S 

( 3p'S (|P) S 

i 

U$) 

But 




- ' 

pW 

p - 


(20) 

where W is the molecular weight and nW = 1- 


\ 2 •' ; 

In addition, 




( H' ) S 

" R^T ( ^S + R^T ( ¥T } S ~ 

pW 

^T2 

(21) 

and 

„VT ‘ T . * 

■ \ -■ 



/ 3W-i 

{ dT } S 

■; <l?)p * <f>T (H>s 


(22) 

•> t; 6 v 1 V ’ 

' w ■ ' ' ' • V.!'. u 

; T 15 tnc 


coefficient 

Then , if we 

; ? * n e d in- ? q c a tier ( 1 S ; . 

define 

' ; 


-• 

°p = ( HV 1> t 


(23) , 

4 


,* - S p 

{ l-D p ) Cp -<l+ D T) 2 Ro 

f 

-# 

(24) 
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Then if r is defined as • • " • ■' : --- 


. - r 



?j> r ° 

(l- D p) C|) -(UDT>2Ro 

f t>e \iozz < $. jf! is vcViS > . 

•r •'•(£) -- • ■ : 

p 



( 25 ) 


(26) 


If we transform Cp into specific heat perunit mass c pm 


then ! 


F = - ' ' : ’ ’ Cpm ' : ' : - 

. : 5 .- . c Pmt 1 - D p>-( 1 + °T) 2 % 




Zd. 


c Pm i 


a. 




(27) 


0 ' 


(28) 


a T 


c Pm 

R^fJ 


1 + Dj ■ 


( 29 ) 


where a.. is the mass fraction of species i and Dy is the 
coefficient defined in Equation (16). 
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From the preceding we conclude that the three thermodynamic 
parameters which must be duplicated are T, ay and cp or 
equivalently Dp, Dy and Cp. 

Unfortunately, the gases entering the nozzle of an actual 
vehicle undergo a substantial change in chemical composition 
as they flow through the duct and this change has. a pronounced 
effect on the values of Dy. and Dp. Since any cold gas simu- 
lation wilT, of neccessity, preclude these chemical reactions, 
it does not appear possible to obtain perfect simulation over 
the entire range of temperatures existing in the nozzle. How- 
ever, gases are available which permit close simulation of r, 

Cp and ay for some conditions of interest. Some of these gases 
are investigated and defined in the next section where their 
variation of critical thermodynamic properties with temperature 
ratio is presented within the temperature range of interest. 



IV. 


COLD FLOW SIMULATION 


Tn order^Lo- as sessj the -r-e-l-atji ve' meri t s-r-of-- var-iou5~g^ses--with 

- I 5 I 1 * j 

regard to their ability to simulate the hot exhaust flow ; 

V ' \ \ • i ■ ; i ! 

through a hypersonic nozzle, fit is first necessary to define 
the properties of a typical hot flow. To do this, we select 

i 

the case of a scramjet propelled vehicle flying at a Mach num- 
ber of 8 >at an altitude of 100,000 feet. The sciramjet engine 
operation is characterized by an inlet static pressure ratio 
of 140 with a total pressure recovery of 60% and; constant pres- 
sure combustion. For hydrogen-air combustion at equivalence 
ratios of 0.6 and 1.0 these operational assumptions result' in 
the following set of gas stream characteristics at the nozzle 
entrance. 


_4_ 

PlP- a .I 

h ( j / kg ) 

p(kg/m 3 ) • 

-t(-°k) 

0.6 

1.56xl0 5 

6x10 s 

0.162 

2821 

1.0 

1.56xl0 5 

4 : . 2 x 1 0 6 

0.204 

2381 


Assuming equilibrium chemistry, the temperature-pressure history 

$ i \ ■ . . • 

through an isentropic expansion from the nozzle ^entrance con- 

ditions to a static pressure' 200 times\lower than the ini-t-jial 

*. ; ■ • \ • ’ | 

value is obtained from the tables of Reference (T) and is {shown 

in Figure (1). It should be noted that the actual, degree of ex- 
pansion would be expected to be s ign i fica-n-t Ty- 1 es s-t-ha-n a •factor 



3000 



FIGURE 1. PRESSURE-TEMPERATURE PATH 'THROUGH EQUILIBRIUM ISENTROPIC EXPANSION. 

■ • • . vr 
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of 200 but th/is wider r,a n g e provides an adequate margjin for 


the purposes ;of co'mpari 


the expansion~curves ofj 

‘I '| 

having also been obtain: 


son 


The 


Figure”! 
ed frorPRe 



ing to 


are pi otted in Figifre (2), 
ife»^ejJ4e 


The values Sioff FT or s'evpr aT“Ry3r bicarb" on gastes “(*0 H 4 
C 2 Hg* C4H3- 2 )|- are shownj in F 1 g ti r e^(3 ') ’ J a h's~of j, tempera-**' 

j t j. Ji 

ture. Since jit is necessary to mjatch the Function r ( ft/ T -j ) 

J j ! ■ | 

I | - | f. , W , | ,? .. ~ — 

through the expansion, Several vaflues of T-j[ were tried for each 



gas and the resulting r* compared with the requirement 


(2). Theses 


of Figure 
the use 


t simulation was found to be obtained by 

ii j- j. j' r 

of C H 4 w i i ft' " i'nXti a 1 Pfemp e ira^^ " o76"”w ^"To 

i * ■ | 

be simulated and by the use of C2H4 with an| initial temperature 
of 500°K when 


Figure ( 4 )v^t ; h|aF th~e ab ijTTty to” provide the proper value of r 


4>=1.0 was, to be simulated. lit is evident from 


is easily attained 




In addi ti on 4 to r, it is also desirable to mptch the variation 

of c p with temperature rati 0 (T/Tk). This variation is pre- 

j j j- • j. 

sented in Figure (5) for the equilibrium f^-air mixtures at 
d> = 0 . 6 and fy=%.0 along wi th that of the meth ine and et iylene d i s j ; 
cussed above. It can be seen fro|i the figure that although thelj 
level of Cp c»n be matched f ai rlyl closely fir 41 = 0.6 s imul ati on , 


, . , . .ftiA—s— -T ~T — r-r-T rrr 

the high initial temperature associated with stoichiometric 


.’»Jr ' 


combustion produces a significant degr&e,.of di ssoci atei on whose 




I: 


effect cannot be closely matched with a relatively cool gas. 


?uiu.qr 


MHvtfnaN-- or 6=011.1 tcnsH-imr r. ihiiojuf NO'mt". .ik 
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FIGURE 4. COMPARISON OF T VARIATION FOR ACTUAL NOZZLE 
GAS AND TWO COOL HYDROCARBONS. 
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FIGURE 5. COMPARISON OF c. p VARIATION FOR ACTUAL NOZZLE GAS 

d ' . r 

• r AND tWO COOL* HYDROCARBONS. (REFERENCE 2) 
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f^inalrly, the value oyf ay- i s ?z.ero -for the nonreactive substi- 
tute gases while it varies from zero to a small negative number 
for the high temperature gases. Although a quantitative evalu- 
ation of the small but finite value of ay for the real gas can- 
no c t be ^mape . he re, it is of >.i nter-est^-to. note, that; dts maximum- -- 
magni tude for <j>=p. 6 is an order of magnitude smaller then it r 
is f , ;fO t r. (j)=l .0;.. ; . v ; . . , i " 


C 


Summarizing, excellant duplication of r can be obtained using 
a cool gas f or L any r vl aue of <{>. Good Cp and ay matching can be 
obtained for moderate values of <j) , but a quantitative defini- 


tion of the error associated with 

w Iv# re .S i's • trve • eagle § i «• 


small mismatches of the three 

„■ a c r t i rt at to r -f f <snh;'..t he 

thermodynamic functions is certainly required. To do this, we 

y s --and wH • is the .local- sontc .velocity ■, gl vee* bj^yEguati.do' 

will integrate the equations of motion used in an equilibrium 

(26T. . from ‘the .definition of. our :-cs.3C<i , fn»fe& ~$y$fc«te;' ,we 

characteristics analysis to analyze a one wave expansion of 


both real gases and substitute hydrocarbons. The results of 
this simple analysis can then be used as a good indicator of 


the error which is associated with the use of a cool substitute 

■( 30 ). . a ad ..( 3 H • o *e r •* ■ < m.i > tep -i n -'ap pt hroae'y:-'- 

gas - ' 

r it . e isvimed .cacstant r ..-ie'ciata tri • 


The momentum equation can be written as 

- ' f 1 ■, „ 2 - ...« 2 «- - - r Z. 

•• % — s~f s Z * ■ V a ’ " C” 


IE. + H dV : 


• V 
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and the isentropic law in the form 


dp _ j, d_p 

i p p 




For a one wave flow, a convenient coordinate system can be de- 
fined by the local Mach line and a normal to it. Defining the 
velocity component along the characterise c as V^., we have for 
one wave flow 


d V t 
dS 


where S is the angle defining the wave orientation from the 
y-axis and a is the local sOTvie/ Velocity, given by Equation 
(26). From the definition of our coordinate system, we have 


V. + a2 


Integrating (30) and (.31) over a small step in Ap through which 
r may be assumed constant, we obtain 


; ; '• t a n ■ n e a • 


(^— 4") a 2 .+ V. + 2 = c 2 

a r ••> § r s i,',*ae d * clr- l j -js erica 1^ s o 1 u t ' •) ■? u *: , .. 


i 1. 1 , P f s ) 2 n d i { ~ , p . i '■ i e ■; c r . b e c’ 


■ n - 


Sub^s trtuti ng :i.nto ( 32 ). r } we ; ha^e ,, 


V T 


. v r u k. ** u * h 


: -uur o . 


■ c r? r o # i r .» * 


C 


(35) 




which yields 



e 2 = e l ~*Y ITT (k 2“ k l } 

(37) 

where 

k s tan~* V t/c 

/l-(V t /c) 2 

(38) 

and 


8 = 0 + p = tt/2 - <5 

(39) 


Equations (33), (34), (37) and (38) can be used to evaluate 
the expansion of a gas through a series of steps in Ap, the 
value of r being a constant in each step. This analysis was 
programmed for numerical solution using the curvefits for 
h(T,p,4>) and r(T,p,<J>) described in Volume I of this report. 

The program listing is presented in the Appendix. The program 
was used to compute the variation of pressure and impulse 
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f uncti on*- wi th flow turni ng- angle through a - one lw ? ave expansion 
for the <j> = 0.6 and 1.0 nozzle entrance conditions described 
earlier. The prbgram was also used fn'th the Ceu rvef i ts of Ref- 

. f. **■* * i I i 


erencehH')-*f;oT'"'tiTe‘ f tmcttf-o h-s^-hT^7'~a-n-d“- f c ffor the gases 

* | i • | « K ’ |i . ' • - * ■’ . ■ 

CH 4 and C 2 H 4 t. Since the temperature: range foiff the hydrocarbons 
was fairly low, the equilibrium isentropic exponent was defined 
by the express i on 


cp/R 

( Cp/R) -1 


(40) 


The results of these calculations are presented in Figures ( 6 ) 
and (7) where it is. clear that fairly good representation of 
pressure distribution and impulse can be achieved with the use 
of cool substitute gases. 




♦Impulse f uncti on., i s .de’fi ned here as', 


<gt 


F/F>i = 


r‘> p 

+ iV 




m 4(r+l) [l + (^-)m 2 ] 


following the notation of Reference (3). 
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■ : • APPENDIX 

PROGRAM LISTING IN BASIC LANGUAGE 


■/ > 




5 : PRINT, '",1 ".."P'S "T'S "RHO 
6 ; -..'.PRINT "GAMMA'S "MACH'S " 

'print 

14 dim ;jcidd l 6'3>Kc iooo] 

1 : 5 : ■' DIM ’BCidOdIVZC 10003 

2D. 7. COTCj 512 

30 ' LET' T=;TQ=1500 ’ ' ' ' 


" > "H '* 
THETA" 


v 


40..: GO.SUB 180 
50 LET E0= (H-Hl >/H 
60 LET HO=H 1 

id :;, let ‘ti=to*ivi • 
oO let ,t-;ti , 

Mi GO.S OB ,180 * * 

Ipo LET 1 El = CH-H 1 > /H 


i'i'o' 

IF ABS( El )< .0001 

THEN 390 

t2o 

LET T9 = TO- EO* C T i 

- TO ) / ( E 1 - EO ) 

T30 

LET E0=E1 


l'40 

LET T0=T1 


ISO 

LET HO=H 1 


160 

LET T1-T9 

'• 

170 

GOTO. 80 :. 

• ' ^ 

180 

LET A=1 . E-07+ ( - • 

1042*Ft2+.8242*F+.987) 


. 190 LET B=.001*< .01 1 67*FT2+ r 1 503*F+.938> 
200 L F.T , C= - . 0 28 4* Ft 2+ . 67 3 1 * F+ i 429 3 
210 IF F < = 1 THEN 250 * ' ' ■ ' '• 




220 LET A=l.E-0 7*Cl'. 78 7*Ft2-5.48*F+5.4> 

230 LET B= » 00 1 * < - • 18 67+ Ft 2+1 . 1 1*F+. 176) 

240 LET C=-.0933*Ft2+3«975*F-2.808 • 

250 I F T,:<= 2000 THEN 330 , 

260 LET A=.000001*< 1 , 792*Ft 2+ . 398 3*F+ . 3 1 ) 

270 LET B=.001*C-9;05*Ft 2-.079l7*F+.24S> 

280 ' LET C=10.86*Ft2-. 1 183*F+.97 
290 IFF < = 1 THEN 330 

300 LET A=.000001*(4.81*Ft2-13.9*F+l 1*59) 

310 LET B=.001*<-23-03*Ff 2+66.82+F-52. 61 ) 

320 LET C=2 7 « 0 5*Ft 2“ 73* 73*F+ 58 • 39 
330 LET Hl=A*Tf 2 + B+T+C 
340 IF T <= 2000 THEN 370 
350 LET 29= . 1 25* CLOC< P) /2 . 3- 5) t 2- . 2 75* (LO G( P) /2 . 3- 5) 
360 LET H1=H1*( 1 + U +F)* < T/2000- 1 ) *29 ) 


370 LET H1=H1*1 .E+06 


380 RETURN 


390 LET M=0 

400 IF T <= 1000 THEN 420 
410 LET M- 2. 1 5E-08* Tt 2+ • 00009 1 * T- • 0 69 5 
420 .LET N= 4 • E-09* Tt 2- « 00002* T- • 0 1 9 
430 IF F <= T THEN 450 ' ' ' 

440 LET, N=. 0 339* SQR < T ) - . 00039 1 * T- . 68 1 
,450 LET G.=^i V833E-07*Tt2 + .pp0p,75*.T+l > 367 
460 ’ IF- T< 500 THEN ' 500. ' j 7 , t' Z * .! . 

470 LET G= 2 • E- 08*Tt 2- • 0001 38* T+l • 423 ! ' 

480 IF T<2000 THEN 500"' : ' J 

490 LET 0=7.267E-08*Tt2-.000457*T+l -85 
•500 LET 0=G+M*CL0GCP)/2.3-5)+N*(F-l ). . 

510 RETURN 7 . ' ' 

512 READ PjF9iRLDLHiViF 

514 DATA 1 . 56E+06> 1 5600> 1 . 62., 5000., 6. E+0 6* 2075> 1 
516 GOTO 600" 

520 PRINT "INITIAL PRESSURE";.', 

530 INPUT P . 

540 PRINT "FINAL PRESSURE"; 

550 INPUT P9 



560 PK1NT "liiij ' ‘ i 

,570 INPUT R 1 7 

580. PRINT "DELTA PRESSURE"; 

590 INPUT. Dl 
600 LET 11=1000 
610 LET JC 13=P 
620 LET KC 1 ] = R 1 
622 GOTO 690 . 

630 PRINT "ENTHALPY"; 

640 INPUT H 

650 PRINT "VELOCITY"; 

660 INPUT V 

670 PRINT "FUEL TO AIR RATIO";" 

680 INPUT F 

690 ' LET H2=H+. 5*V+V 

700 GOSUB 30 

710 LET C8=P/RltG 

,720 LET R=R1 ; > J 

722 GOSUB 2000 

724 LET ZC 1 3=ATN< Vl/Cl/SQRC 1- C VI /Cl > t 2> > 
726 LET BC 1 3=ATN< 1/SQRC < V/Al > t 2- 1 ) > 

728 LET T1 = 0 
730 FOR 1 = 2 TO I 1 
740 LET C8 = P/Rt 0 
744 IF I <2 THEN 750 
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746 : BETAC0.«) = #£##j?#. MU=#««# THETA 

750 IF INT(I/20)*20 <> I THEN 770 - 

760 PRINT I,P/ T/R/H . 

'765' PRINT G/ V/Al, 57. 3*T1 
766 PRINT 

770 let jci3=jci-13-di 

780 LET KCI3=<JCI3/C8>t(l/G> 

'790 LET V=SQR<V*V+2*G/(G-i:)*CJCI-13/KCI-13-JCI3/K'CI3)> 
800 LET H=H2~ * 5* V*V 

810 LET P=JC II * 

820 LET R=KC II 

830 GOSUB 30 ' . 

835 GOSUB 950 i 

840 LET 12=1 

850 IF JCI 3 <= P9 THEN 890 
860 IF JC I 3 -D1>P9 THEN 880 
870 LET Dl = JC I 3-P9 

880 NEXT I / 

890 LET I=KC I23*V*V+JC 123 ' ‘ 

900 - LET R=KC 1 23 , • ' - 

910 LET P=JC 123 

920 GOSUB 30 ' _ 

925 GOSUB 950 

930 PRINT I2>JCI23jTjKCI23>Hj G> Vj I ' 

940 STOP , - ' 

950 GOSUB 2000 

990 LET ZCI3=ATNCV1/C1/SQRC 1-CV1/C1 )«2>) 

1000 LET BCI3=BCI-13-SQRCG1)*CZCI3-ZCI-13) 

1010 LET U=ATN( 1 /SQRC ( V/Al.) t 2- 1)3 ~ 

1020 LET T1 = BC I 3 - U 
4024 GOTO 1030 

1025 PRINT USING 746/ I / BC I 3 / 57. 3*U/ 57- 3*T1 
1030 RETURN 
1990 STOP 

2000 LET A1=SQR< C*P/R> 

2010 LET G1=<G+1 l/CG-l > ' 

2020 LET C1=SQR( Gl+Al t 2+Vt 2-A1 t 2) 

2030 LET V1 = SQR(V»2-Alt2) 

2040 RETURN 

3000 END - ' 



